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Abstract 

It is shown that the perturbative expansions of the correlation 
functions of a relativistic quantum field theory at finite temperature 
are uniquely determined by the equations of motion and standard 
axiomatic requirements, including the KMS condition. An explicit 
expression as a sum over generalized Feynman graphs is derived. The 
canonical formalism is not used, and the derivation proceeds from 
the beginning in the thermodynamic limit. No doubling of fields is 
invoked. An unsolved problem concerning existence of these pertur- 
bative expressions is pointed out. 
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1 INTRODUCTION 



The traditional way of describing thermal equilibrium states of an infinitely 
extended quantum system, in particular of a quantum field theory, begins by 
restricting the system to a finite volume V, defining the canonical or grand 
canonical equilibrium by means of the familiar density matrices, and then 
going to the limit V — > oo (the "thermodynamic limit") for the quantities 
for which this limit can be expected to exist [1,2]. This applies especially 
to the correlation functions of the fields and closely related objects like the 
expectation values of time ordered field products. Up to now most actual 
calculations of such functions have been based on this approach, using a 
Hamiltonian or Lagrangian formalism at finite V. 

Another description of equilibria and their local disturbances, which can 
be used directly in the thermodynamic limit, has been developed in the frame- 
work of the algebras of local observables [3]. In this approach equilibrium 
states are characterized through an analyticity requirement for correlation 
functions, the so-called KMS condition. In the present paper we intend to 
show that this axiomatic method, suitably adapted to a field theoretical con- 
text, is perfectly capable of handling dynamical problems. More exactly, it 
will be shown that perturbative expansions for the correlation functions of 
a relativisitic field theory, and related functions, can be derived directly in 
the thermodynamic limit, not making use of the canonical formalism, but 
using as only inputs the equations of motion and the axiomatic requirements 
that the correlation functions must satisfy. The result is represented as a 
sum over generalized Feynman graphs. For the special case of time ordered 
functions it agrees with the well-known result of the canonical approach. 

Dispensing with the canonical formalism is also a major difference be- 
tween our approach and thermo field dynamics [2,4], a Fock space method 
developed by H. Umezawa and coworkers. We differ also from thermo field 
dynamics by not invoking a doubling of fields, and by not assigning a basic 
role to particles, including quasi-particles. We hold particles to be secondary 
objects of the theory, of great phenomenological importance, but little funda- 
mental significance. In this respect we differ also from the views put forward 
by Landsman in ref. [5]. 

We consider only the ^-model, i.e. the theory of a scalar hermitian field 
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<3>(x) satisfying the equation of motion 



(□ + m 2 )$ 
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iV($ 3 ), 



(1.1) 



where N stands for a normal-product prescription taking care of renormal- 
ization. The restriction to $4 is merely a matter of convenience. The gen- 
eralization of the method to other models, including gauge theories in local 
gauges, is straightforward. 

We are interested in the correlation functions 



where (•) denotes the expectation value in a thermal equilibrium state with 
temperature T > 0. These correlation functions describe the full physical 
content of the theory: all observable quantities can in principle be derived 
from them (for examples see e.g. [1,2]). This is so because knowledge of the 
W allows the reconstruction of the full representation of the field algebra by 
means of the GNS construction [2,3], which yields a Hilbert space represen- 
tation of the field algebra with a cyclic vector |), such that (A) = (\R(A)\) 
for any sufficiently regular function A($), with R(A) its representative. All 
observables of the theory are supposed to be of this form, and local distur- 
bances of the equilibrium are created by applying suitable functions F(&) to 



More generally, we consider the set of functions (or rather, distributions) 



Here the X a are non-overlapping sets of 4- vectors Xi, the s a are signs, and 
T ± (X) denotes respectively the time-ordered or anti-time-ordered product of 
the fields $(x;), Xi E X. If each X a contains only one variable, then W is the 
correlation function W(x±, . . . , x^) irrespective of the choice of the signs s a . 
For Af=lwe obtain the usual time-ordered and anti-time-ordered functions 
(Green's functions) of the theory. In the sequel the signs s a will be frequently 
suppressed when they are not essential to understanding. 

In a previous work [6], henceforth quoted as V, we have derived perturba- 
tive expressions for the functions W in terms of generalized Feynman graphs 
for the case T = 0, in which case (•) denotes the vacuum expectation value. 
This derivation uses neither the Hamiltonian nor the Lagrangian formalism, 



W(x!,...,x n ) = ($(xi) ...$(x n )), 



(1.2) 
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W(X U Sl \... \X N , s N ) = (T S ^X 1 ) . . . T SN (X N )). 



(1.3) 
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but is instead relying on the Wightman axioms as an essential input. We 
propose to generalize this method to the case T > 0. Our way of proceeding 
is closely modeled on that taken in V. The ideas and results of that paper 
will be freely used. Equation (n.m.) of V will be referred to as eq. (V.n.m.). 
The method conists in solving the differential equations for W which follow 
from the field equation (1.1), by a power series expansion in the coupling 
constant g, using the axiomatic properties of W as subsidiary conditions. 

The paper is organized as follows. The assumptions on which the for- 
malism is based will be stated and discussed in section 2. An explicit formal 
expression for the W's as sums over generalized Feynman graphs will be 
stated in section 3 and shown to possess the required properties. In section 
4 the ultraviolet (UV) divergences of these graphs will be removed by renor- 
malization. It will be pointed out, however, that renormalization does not 
guarantee the existence of the resulting finite-order expressions, on account 
of the local singularities of the integrands. This problem remains unsolved. 
Finally we will show in section 5 that the expressions of sections 3 and 4, 
provided they exist, are the only ones satisfying the assumptions stated in 
section 2. 

2 ASSUMPTIONS 

In this section the assumptions on which our derivations are based, in addi- 
tion to the field equation (1.1), will be enumerated. 

The following conditions for W and W are taken over unchanged from V. 

a) The W and W are invariant under space-time translations and under 
space rotations. Invariance under Lorentz boosts cannot be demanded for 
T > 0.. 

b) Locality holds, i.e. W(X) is invariant under the exchange of two neigh- 
bouring variables Xi,x i+i , if x^ — x i+ i is space-like. 

c) The reality condition 

W(X U Sl \... \X N , s N )* = W(X N , -s N \ . . . \X U - Sl ) (2.1) 

holds. 

d) The functions W are permutation invariant within each sector X a , and 
they satisfy the splitting property 

W(. . . \X ± U X 2 , +| . . .) = W(. . . \X U +\X 2 , +| . . .) (2.2) 
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if x® > for all Sj 6 Xi, Xj G X2. This property shall hold in every Lorentz 
frame, not only in the rest frame of the infinite system under considera- 
tion. To that extent we retain Lorentz invariance. These conditions are not 
merely a definition of time ordering. That they can be satisfied is intimately 
connected with locality. 

The cluster property needs a more careful discussion than it was accorded 
in V. Let X, Y, be two non-overlapping sets of 4- vectors. Then the cluster 
property states that 

lim W(X, Y + a) = W(X)W(Y) , (2.3) 

|a|— +00 

where Y + a means that all vectors in Y are translated by a, and a tends to 
infinity in a space-like direction. Since derivation with respect to the coupling 
constant need not commute with the limit in (2.3), we cannot expect eq. (2.3) 
to hold separately in each order of perturbation theory, except in the lowest, 
free, order. But perturbation theory can also be viewed as an expansion in 
powers of H. And we can demand that eq. (2.3) hold for each W(X) in the 
lowest nonvanishing order in k. This is all that will be needed to establish 
uniqueness. We set c = 1, so that time and space have the same dimension. 
The dimension of the field $ is then [mass/length] 1 ^ 2 . In order to make the 
field equation (1.1) dimensionally consistent, the coupling constant in the 
interaction term reads g/h, with a dimensionless g. We shall nevertheless 
set H — 1 in our equations, and merely note the correct ^-exponents at the 
points where they are essential. 

The spectral condition of the vacuum representation does not hold at 
positive temperatures. It is replaced by the KMS condition, which we use 
in its p-space form (see [3], lemma 1.1.1 of chapter V). Define the Fourier 
transform $(p) of $ by 

$(p) = (2n)- 5/2 J dxe tpx $(x), (2.4) 

and the Fourier transform VV(Pi| . . . |Pjv) of W(Xi| . . . \X^) accordingly. Let 
Pi, ... , P/v, Qn+i, ■ ■ ■ , Qm, M > N, be finite, non-empty sets of 4- momenta. 
Define 

P°= E Vl (2.5) 
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Then the KMS condition states that 

W(Pi\ ■ ■ • I^atIQat+iI . . . \Qm) = eP p °W(Q N+1 \ . . . \Qm\Pi\ ■ ■ ■ \P N ), (2.6) 

where f3 — is the inverse temperature. 

The normalization conditions stipulated in V are taken over unchanged 
for the case T = 0. They are the standard conditions demanding that m be 
the physical mass of the particles associated with the field $ at T = 0, and 
defining the coupling constant g and the field normalization in terms of the 
Green's functions. In addition we demand that <£>, considered as an element 
of an abstract, representation independent, field algebra, is T-independent. 
This implies that the field equation is T-independent, meaning that the pa- 
rameters m and g and the subtraction prescription N do not depend on the 
temperature. The prescription N is assumed to be the conventional one. In 
particular it should not contain oversubtractions which would destroy the 
renormalizability of the theory. In perturbation theory the T-independence 
of N means that the UV subtractions and the subsequent finite renormal- 
izations have the values used at T = for all T. Note that in the BPHZ 
method [7], which we will be using, the subtractions do not involve prior reg- 
ularizations or any formal juggling with divergent quantities. The procedure 
is therefore well defined. For more details we refer to section 4. A nonper- 
turbative definition of N can possibly be given with the help of a Wilson 
expansion, as explained in section IV. 2. of ref. [7], i.e. by a point-splitting 
method with suitable singular coefficient functions. T-independence of N 
could then be defined as T-independence of these coefficients. We will not 
explore this possibility any further. 

Note that the parameter m is for T > not the physical mass of a particle. 
Quasi-particle masses are determined by the positions of the singularities of 
the clothed propagator, which are T-dependent. 

3 UNRENORMALIZED SOLUTION 

The coefficient W a (X) of order a in the power series expansion 

oo 

W{X) = Y J 9 a W a {X) (3-1) 

cr=0 
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is determined as a solution of the system of differential equations 

(□* + m 2 )(. . . ...)«. = —<••■ iV(<f fe) 3 ) • • •).-! (3-2) 

satisfying the subsidiary conditions described in section 2. For <r = the 
right-hand side is zero, for o > it can be calculated if the problem has 
already been solved in order a — 1 . 

We first state a formal, unrenormalized, graphical representation for 
W (T (Xi| . . . \X N ) and then show that it does indeed satisfy all the require- 
ments. 

For the free propagators we use the following notations. We define first 
the p-space expressions 

1 i 1 

(27r) d (2tt) p — ^ + 2e 

which are the usual vacuum propagators, and their thermal extensions 

D + (p) = A + (p) + C(p), D F (p) = A F (p) + C(p) (3.4) 
with the T-dependent correction term 

^)=(^^^ 2 "- 2 )- (3-5) 

Here 5+(p) = 9(p )S(p 2 — m 2 ) is the 5-function of the positive mass shell, and 
u; = (p 2 + m 2 ) 1 / 2 . Note that the additional term C is the same for D + and 
D F . 

The a;-space versions of these propagators are defined by 

A. (x)=if d 4 P A. {p)e~ ipx , D. (x) = i J d 4 P D. {p) e - ipx . (3.6) 

The thermal propagators (3.4-5) need not be taken over from the tradi- 
tional formalism. They can be derived in our framework by the methods 
that will be used in section 5 to establish uniqueness of our solution. But in 
the present section we rely on an a posteriori justification, by showing that 
these forms give rise to expressions with all the required properties. 

W CT is represented as a sum over generalized Feynman graphs which are 
defined as follows. Draw first an ordinary Feynman graph of the <3> 4 theory 
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with \X\ — J2 a \X a \ external and a internal vertices. Here \X a \ is the num- 
ber of points in the set X a . The graph need not be connected, but must 
not contain any components without external points. This graph is called 
the "scaffolding" of the generalized graph. Next, it is partitioned into non- 
overlapping subgraphs, called "sectors", such that the external points of a 
set X a belong all to the same sector, but variables in different X a to differ- 
ent sectors. There may also exist "internal sectors" not containing external 
points. The sectors are either of type T + or T~. For external sectors this 
sign is given by s a . To each sector S we affix its number u(S) according to 
the following rules. 

i) v(S) = a for the external X Q -sector. 

ii) If s a = s a+ i there may be an internal sector with number a + |. Its 
type is the reverse of the adjacent external sectors: s a+ i = —s a . If 
s a = — s a+ i no such intermediate internal sector exists. 

iii) If si — sn there may be an internal sector with number N + | and 
s N+ i = — sjv- Equivalently we could give this additional sector the 
number |, but we must choose one of these possibilities and use it 
consistently. 

With a partitioned graph we associate a Feynman integrand, which we 
state first in x-space. To the external points correspond the external vari- 
ables X; t . To each internal vertex we assign an integration variable Uj, j = 
I,..., a. Zi denotes a variable which may be either external or internal. 
Within a T + sector the usual Feynman rules hold: each internal vertex car- 
ries a factor —ig, a line connecting the points Zj and Zj carries the propagator 
—iDp(zi — Zj). In a T~ sector the complex-conjugates of these rules apply. 
A line connecting points z iy Zj in different sectors, with Zi lying in the lower- 
numbered sector, carries the propagator —iD + (zi — Zj). The graph must 
be divided by the usual symmetry number if it is invariant under certain 
permutations of points and lines. 

In p-space, integration variables are assigned to the lines. The vertex 
factors in factors are ^f(2n) 4 g and the propagators are D F (p) or (D F (p))* 
respectively. Lines connecting different sectors carry propagators D + (p), the 
momentum p flowing from the lower to the higher sector. For each external 
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point there is a factor (27r) 3 / 2 . Momentum is conserved at each internal 
vertex. 

For T = the rule ii) governing internal sectors seems to differ from the 
corresponding rule in V, where instead of one (possibly empty) intermediate 
sector of different type we had chains of intermediate sectors of the same type 
as the bracketing external sectors. However, it can be shown with the help of 
Lemma V.3.1 that the two formulations are equivalent. Graphs containing 
non-empty internal sectors according to iii) vanish for T = 0, wherefore they 
do not occur in V. 

For the Green's functions (N = 1) our rules agree with those of the con- 
ventional real-time formalism. Our graphs are equal to those of the Keldysh 
formulation [8]. 

The first point to be checked is that the above prescription gives an 
unambiguous result for the correlation functions W. The problem is that 
in the definition W(xi, . . . ,x n ) = W(xi, Si| . . . \x n , s n ) the signs s a can be 
chosen arbitrarily. But it can be shown that the sum over all graphs with 
the same scaffolding does not depend on the choice of these signs. The proof 
is modeled closely on the corresponding proof in V and will only be briefly 
indicated. Consider a given scaffolding. Let S denote an internal T + sector, 
S x an external T + sector with external variable x, S and S x T~ sectors. 
A product S X S or the like denotes the sum over all partitions of a given 
subgraph into two adjacent sectors S x , S. The propagators connecting these 
two sectors are included. We will first show that the choice of sn is irrelevant. 
The graphs belonging tos7v = + orsAr = — respectively differ only in the 
sectors with number u(S) > N — 1. In the case sat_i — +,Si — — these 
variable sectors are (we put x N = x) S x + SS X for = +, and S x + S X S for 
sn = —■ But 

S X + SS X = S x + S X S 

is a special case of Lemma V.3.1, which remains valid for our new propaga- 
tors. Similarly, if s± = sjv-i = + we must show that 

S X + S X S + SS X + SS X S = S x , 

or 

(S x + S X S -S x - SS X ) + S(S X + S X S -S x - SS X ) + (S + S + SS)S X = 0, 
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which is correct because all three brackets vanish as a result of Lemma V.3.1. 
The proof of independence on the other s a follows similar lines. The case 
1 < a < N is considerably simplified compared to V by the new formulation 
of the rules concerning internal sectors. 

Of the conditions stated in section 2 invariance under translations and 
rotations, and the symmetry of W within T ± factors are trivially satisfied. 
That the equation of motion (3.2) is satisfied, is shown exactly as in V. At 
the present formal level the product iV($ 3 ) denotes Wick ordering : $ 3 (a;) : 
with respect to the free vacuum field. 

The proof of the reality condition (2.1) and Ostendorf's proof [9] of the 
splitting relation (2.2) can also been taken over unchanged from V. Since 
the non-invariant part C(x — y) is the same in D + and in Dp, the crucial 
relation D F (x — y) = D + (x — y) if x° > y° holds in every orthochronous 
Lorentz frame. Hence the same is true for the splitting relation. From this 
and the symmetry within sectors we can prove locality. Let (x — y) 2 < 0. 
Then there exist two Lorentz frames such that in one x° > y°, in the other 
x° < y°. Hence 

W„(...,x,y,...) =W ff (...\x,y,+\...) = 
= W a (...\y,x,+\...) = W a (...,y,x,...). 

Our propagators, being the 2-point functions of a free field, contain for 
dimensional reasons a factor H (this refers to x-space). The vertices carry a 
factor ft -1 , so that the 2n-point function of order a contains a factor K a+a . 
The term of lowest order in h is therefore the free a — expression. It 
is simply a sum over all possible products of free 2-point functions, and 
obviously satisfies the cluster property, because D + (£) and Dp(£) converge 
to zero if £ tends to infinity in a space-like direction. Hence our weak version 
of the cluster property is satisfied. 

Yet to be proved remains the KMS condition (2.6). It is easy to see 
that a graph contributing to the left-hand side of (2.6) becomes a graph 
contributing to the right-hand side, if the directions of all lines connecting 
sectors with v(S) < N+l to sectors with u(S) > N + l are reversed, and vice 
versa. This means that the corresponding propagators D+(k) are replaced by 
D + (—k). From the definitions (3.3-5) we find D+(k) = e /3k °D + (—k). Hence 
the two variants of the considered graph differ by the factor exp(/3J2 a )> 
the sum extending over the momenta of the lines in question. But K = Y,k a 
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is the total momentum flowing from the P-part of the graph to the Q-part, 
so that K° = P° as defined in (2.5). As a side remark we note that a similar 
argument can be used to show that it is immaterial whether extremal internal 
sectors are assigned the number \ or TV + \. 

4 RENORMALIZATION AND THE EXIS- 
TENCE PROBLEM 

As yet, the expressions derived in the preceding section have only a formal 
meaning. There remains the question of the existence of the integrals sym- 
bolized by the graphs. 

The UV problem is concerned with the behaviour of the integrands at 
infinity in momentum space. It can be handled exactly as was done in V 
for the vacuum case. We note first that loop integrals over loops extending 
over more than one sector are finite because of the strong decrease of D + (k) 
for k — > — oo and momentum conservation: primitively divergent subgraphs 
exist only within sectors, where they can be treated with conventional meth- 
ods. We choose BPHZ renormalization (see [7]), which introduces suitable 
subtractions in the integrand, before integrating, thereby avoiding the need 
for regularization. The subtractions are found by expanding the integrands 
of potentially dangerous subgraphs, the "renormalization parts", into power 
series of sufficiently high degree in their external momenta. The tempera- 
ture dependent part C of the propagator Dp decreases exponentially fast 
at large momenta and is innocent of any UV problems. We can therefore 
define the mentioned subtraction terms using Ap instead of Dp inside the 
renormalization parts, without destroying the UV convergence achieved by 
the subtraction. We also define the finite renormalizations of the BPHZ pre- 
scription to be T-independent, giving them the values needed to satisfy the 
normalization conditions at T = stated in V. These rules are the expression 
of the required T-independence of the renormalization prescription N in the 
field equation (1.1). An important effect of this limited subtraction is the 
emergence for T > 0, of lines connecting a vertex to itself. But they carry 
the propagator C(k), hence the loop integral over k exists, and is the same 
in T + and in T~ sectors. 

But the UV problem is not the only existence problem we are faced with. 
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There is also the problem of the local singularities of the integrand, again in 
p-space. The integrand is a product of distributions in variables which are 
quadratic functions of a complete set of independent external and internal 
momenta. The variables in a given subset of propagators may be dependent 
(i.e. their gradients in momentum space may be linearly dependent) on cer- 
tain manifolds, in which case the product of propagators does not necessarily 
define a distribution. It is no longer clearly integrable even in the sense of 
distributions. Contrary to assertions found in the literature, the problem 
is not restricted to the case of two propagators separated by a self energy 
insertion, and thus both depending on the same variable. As an example of 
a more complex situation, consider the integral / dkD.(k)D ■ (p — k) over a 
two-line loop, the dots standing for either + or F. The mass shell 5's in the 
two factors coalesce at p — 0, hence the integral diverges at that point, and 
this singularity in the external variable p (external to the considered 2-line 
subgraph) is not removed by renormalization. Closer inspection shows that 
the singularity is of first order. A chain of n two-line bubbles will then pro- 
duce a singularity of order \p\~ n , which is not integrable for n > 4. Nor is it 
defined in another way as a distribution. The remaining integration over p, 
which may be an internal or external variable of the full graph, is therefore 
not defined: individual graphs containing such chains diverge. These diver- 
gences may cancel between graphs with the same scaffolding but different 
sector assignments of the vertices in the chain. But we are not aware of a 
proof to this effect, even for this still relatively simple example. 

How is this problem solved in the vacuum case? This is easiest for the fully 
time ordered function f(P) = W(P, +). If m > this function is everywhere 
the boundary value of an analytic function in complexified variables pi. The 
same is true for the Feynman propagator (k 2 — m 2 + ie)~ l . For variables 
{pi} in the domain of analyticity of f we can deform the integration contours 
for the internal variables kj into the complex in such a way that they never 
meet the singularities at k 2 = m 2 . The integrand is then a smooth function, 
and there are no problems with the local existence of the integral. But the 
presence of a <5-term in Dp destroys analyticity, so that the method does 
not work for positive T. For m — 0, Ap(k) is singular at k — and is 
there not a boundary value of an analytic function. The same holds for f at 
points where a partial sum of p^s vanishes. These singularities can lead to 
infrared divergences, which need special attention. An x-space method for 
achieving this, and also for proving existence of the general W a , has been 
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described in V. It relies heavily on the fact that A+(£) is analytic in im£° < 
and decreases at least of second order for im£° — > — oo. But D + (£) is only 
analytic in the strip —(3 < im£° < 0, so that this method can also not be 
extended to positive T. 

At present the problem of existence of W CT remains unsolved. 

5 UNIQUENESS 

It will be shown that the expansion described in sections 3 and 4, assuming 
its existence in finite orders, is the only expression satisfying all the require- 
ments. 

Assume that uniqueness has been established up to order o — 1. Let 
W*,W 2 , be two solutions of the equations (3.2) plus subsidiary conditions. 
Then their difference 

h a ( Xl , ...,x n ) = Wfau ...)- W 2 (x ± , . . .) (5.1) 

satisfies the homogeneous equations 

(D i + m 2 )h a (...,x i ,...) = (5.2) 

and all the subsidiary conditions. We want to derive the most general form 
of h a . 

The Fourier transform h a (pi, . . . ,p n ) must contain a factor (pf — m 2 ) for 
each i. 
Define 

K(. ..[x,y]...) = h a (...,x,y...)- h a (. ..y,x...) (5.3) 

and let h a (. . . [p, ?]...) be its Fourier transform. As in V we can prove that 
h a (. . .[p,q] . . .) has its support contained in the manifold p + q = 0, and must 
therefore be of the form 

h a {. ■ ■ \p, ?]..•) = 5\V + q)S( P 2 - m 2 ){/ + (p) + e(po)/-(p)} , (5.4) 

where the dependence on the variables indicated by dots has been suppressed. 
By locality the support of the function 6(x° — y°)h a (. . . [x, y] . . .) is contained 
in the set (x — y)eV + , the closed forward cone. Its Fourier transform 



Mp) p» +/ _ (p) . i 



^(p) (Po + ^ e ) 2 — P 2 — m2 (Po + ^ e ) 2 — p 2 — rn 
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considered as a function of p + q and q, is in p analytic in the forward 
tube {Imp e V+} and is polynomially bounded in the slightly smaller tube 
{(Imp — a) E V + } for any a E V + [10]. The p -dependent factors in the 
above expression have the correct analyticity. In order for this to be the case 
for the full expression, the p-dependent terms u~ 1 f + and /_ must be entire 
functions of p. Furthermore, they must be polynomially bounded, hence 
they are polynomials. As a result we find 

K{. ..[p,q]...)= 5\p + q)5(p 2 - m 2 )Mp)F( P ) + e(po)G(p)} (5.5) 

with F and G polynomials in p. F and G may also depend on the variables 
not shown explicitly, and factors only depending on those other variables 
have been omitted. 

Consider next the double-commutator function 

h a (. . . [[x,y],z] ...) = h a (. . . [x,y]z.. .) - h a (. ..z[x,y] . . .) ■ (5.6) 

By locality its support is contained in the set S = {(x — z) 2 > or (y — z) 2 > 
0}. But by transforming (5.5) into x-space we find that the right-hand side 
of (5.6) depends on x and y only in the combination x — y, i.e. it is invariant 
under simultaneous translation of x and y by the same 4-vector. But no 
subset of S, apart from the empty set, is invariant under such a simultaneous 
translation. Therefore we find 

h a (...[[x,y]z]...)=0 

Finally we find from the KMS condition (2.6) that 

(1 -exp(-[3p° 1 ))h a (p l ,...,p n ) = h a (p 1 ,...,p n )-h <T (p 2 ,...,p n ,Pi) 

n 

= X)MP2,...[pi,Pi]...,Pn), (5.8) 

i=2 

from which h a can be obtained through division by (1— exp(— Pp®)). Reinsert- 
ing the resulting expression in the terms on the right-hand side of eq. (5.8), 
and iterating this procedure sufficiently many times, we find that h a (p\, . . . , p2n) 
must be of the form 

h a (pi, ■ ■ -,P2n) 

n 

E nPn, , P ij II [1 - eM-frlT'Sipl ~ m 2 )S\p la MPS) 

a=l 



(5.7) 
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where the sum extends over all partitions of the 2n variables into n ordered 
pairs (Pi a ,Pj a ), i a < ja, and F is a polynomial in and cu(pi). The odd- 
point functions vanish in $ 4 theories. 

Now, if F is a genuine polynomial, not a constant, then the expression 
(5.9) shows a bad high-energy behaviour for p® a — > oo, which would destroy 
renormalizability. Hence F must be a constant, and h a has in x-space the 
form 

h a (xi, . . . , x 2n ) = c an Ii D +( x ic ~ x jJ- ( 5 - 10 ) 

pairings a 

For dimensional reasons c an must contain a factor H n , independently of a. 
Insertion of such expressions in the right-hand side of the equation of motion 

(3.2) leads to terms of higher order in h. This means that the homogeneous 
terms (5.10) are the contributions of lowest order to the 2n-point functions. 
Since their functional form does not depend on a, their contributions to 
W(xi, . . . ,x 2n ) can be summed over a, in the sense of formal power series, 
to yield 

c n (<?) £n £>+(•••) (s.ii) 

a 

as term of lowest order in H. The normalization condition for the 2-point 
function determines the value of C\ for T = to be C\ = —i. With a different 
choice of C\ for T > our T-independent definition of the normal product iV 
would no longer remove the ultraviolet divergences, i.e. such a choice would 
destroy the validity of the field equation (1.1). Hence we must have 

ci = -i (5.12) 

identically in g and T. 

By our assumptions the expression (5.11) must satisfy the cluster property 

(2.3) . Using induction with respect to n one finds the unambiguous result 

c n = H) n • (5-13) 

This value does not depend on g, which proves the desired unicity 

h a = (5.14) 

for a > 0. In the lowest order a = 0, ho is the only nonvanishing contribution 
to the 2n-point function. It has the conventional free-field form, as has 
already been anticipated in section 3. 
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It is a curious aspect of this argumentation that it does not work in 
the free case g = 0, where renormalizability is not at stake. For the free 
theory one could arrive at the same results faster by postulating the usual 
canonical commutation relations. For interacting fields this method is rather 
less convincing, because in that case the CCR's have no rigorous meaning, 
since interacting fields cannot be restricted to sharp times. 
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